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ABSTRACT 

We constrain the total mass distribution in the cluster A3571, combining spatially resolved ASCA 
temperature data with ROSAT imaging data with the assumption that the cluster is in hydrostatic 
equilibrium. The total mass within r 50 o (1-7 h^ Q Mpc) is M 50 o = 7-8-2.2 x 10 14 ^50 -^0 a ^ 90% confidence, 
1.1 times smaller than the isothermal estimate. The Navarro, Frenk & White "universal profile" is a 
good description of the dark matter density distribution in A3571. The gas density profile is shallower 
than the dark matter profile, scaling as r~ 2A at large radii, leading to a monotonically increasing gas 

mass fraction with radius. Within r 500 the gas mass fraction reaches a value of f gas = ^-^-am ^50^ 
(90% confidence errors). Assuming that this value of f gas is a lower limit for the the universal value of 
the baryon fraction, we estimate the 90% confidence upper limit of the cosmological matter density to 
be OT < 0.4. 

Subject headings: cosmology: observations - dark matter - galaxies: clusters: individual (A3571) - 
intergalactic medium - X-rays: galaxies 



1. INTRODUCTION 



Measuring cluster masses has significant implications for 
coamology. Aaauming that the cluatcr maaa content rep 



resents that of the LJnivorsc, measured total and bary - 



onic mass distributions in a cluster, combined with the 
Big Bang nucleosynthesis calculations and observed light 
element abundances, can be used to constrain the cosmo- 
logical density parameter (White et al. 1993). A measured 
cluster mass function could be used to constrain cosmolog- 
ical parameters via the Press-Schechter formalism. How- 

■ ever, the mass cannot be observed directly since most of 
the cluster mass is dark matter, and one must rely on ob- 

\ served quantities like gas temperature or galaxy velocities 
and assume the state of the cluster to derive the mass. 

In this paper, we estimate the total mass for the A3571 
cluster under the assumptions of hydrostatic equilibrium 
and thermal pressure support. Until recently, most hydro- 
static X-ray mass estimates have been made assuming that 
the gas is isothermal at the average broad beam tempera- 
ture. However, the total mass within large radii is only as 
accurate as the local temperature at that radius. ASCA 
observations provide spatially resolved temperature data 
for hot clusters and yield their 2D temperature structure. 
A large number of ASCA clusters show s that the temper- 



ature declines with increasing radius (Markevitch et al 



1998), in qualitative accordance with hydrodyna.mic clus- 



ter simulations (e.g. Evrard et al. 199C , Bryan fe Normar] 



1997 



Burns et al. 1999). This implies that the total mass 



within small radii is greater, while at large cluster radii it 
is smaller than that derived assu ming isothcrmality, which 
has also been o bserved in A2256 ( Markevitch fe Vikhlinin| 
1997b| , A20 29 ([Sarazin, Wise fc Markevitch 1998]) , A496 



and A2199 (Markevitch et al. 1999| ) and A401 ( [Ncvalaincr] 
et al. 1999a). Consequently, the gas mass fraction within 



large radius where the cluster may be a fair sample of 
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the universe, is larger than that derived assuming isother- 
malit y which further aggravates the "baryon catastrophe" 
(e.g. IWhite et al. 1993|, |White fc Fabian 1995|, |Ettori & 
Fabian 199H| , |Mohr et al. 1999| ) 



A3571 (z = 0.040) is suitable for measuring the dark 
and total mass distributions, since it is bright and hot (~ 
7 keV) allowing accurate temperature determinations with 
ASCA. In deed, the A3571 tempe rature profile used in this 
work (see Markevitch et al. 1998| ) is among the most accu- 
rate for all hot clusters. The three ASCA pointings cover 
the cluster to rsoo (the radius where the mean interior den- 
sity equals 500 times the critical density, approximately 
the radius inside which hydrostatic equilibr ium holds, ac- 
cording to si mulations of Evrard et al. 1996 ). A3571 has a 
cooling flow (Peres et al. 199? ), but it is weak enough not 
to introduce large uncertainties in the temperature deter- 
mination. 

We use Hq = 50 /150 kms -1 Mpc -1 , Q = 1 and report 
90% confidence intervals throughout the paper. 

2. ROSAT ANALYSIS 

We processed the ROSAT data, consisting of a PSPC 
pointing rp800287 , using Snowdcn's So ft X-Ray Back- 
ground programs ( [Snowden et al. 1994 ), which reduced 
the total exposure by 25% to 4.5 ks. The spatial analy- 
sis was restricted to the energy band of 0.73 - 2.04 keV 
(Snowden's bands R6-R7) to improve sensitivity over the 
X-ray background. The surface brightness contour map 
(smoothed by a Gaussian with a = 1') is shown in Figure 
nk. The data show no obvious substructures and no devia- 
tions from azimuthal symmetry, except for a slight elliptic- 
ity. Fabricant et al. (1984) showed that for A2256, whose 
X-ray brightness distribution is more elliptical than the 
one of A3571, the true elliptical total mass is very close to 
the hydrostatic total mass derived assuming spherical sym- 
metry. Furthermore, Vikhlinin et al. (1999) divided the 
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X-ray total mass estimate for A3571 



PSPC data of A3571 into several sectors, fitted the bright- 
ness with an azimuthally symmetric model and found that 
the azimuthal variation of the gas density gradient (due to 
ellipticity) was 6% of the global value, which would indi- 
cate a similar error in the total mass, which is negligi- 
ble compared to total mass errors obtained with spherical 
model. Therefore, in the following analysis we assume the 
cluster to be azimuthally symmetric. 

We excluded point sources and generated a radial sur- 
face brightness profile in concentric annuli of width rang- 
ing from 15" at the center to 7' at a radial distance of 
40'. In the radial range of ASCA pointings (r < 35'), we 
included only the ROSAT data from the sky areas covered 
by ASCA. A cooling flow with a mass flow rate of 40 - 130 
Mo/yr and a cooling radius of 1.2 - 2.2 arcmin has been 
detected in the center of A3571 by Peres ct al. (1998). 
Our data show a significant central brightness excess over 
the (3- model (Figure ^), in agreement with the reported 
cooling flow. Therefore we excluded the data within the 
central r < 3' from the fit. 

We fitted the observed profile with a (3- model 



1(b) 




(-3/3+i) 



CXRB 



(1) 



( Cavaliere fc Fusco-Femiano 1976 ), where b is the pro- 
jected radius. We fixed the cosmic X-ray background 
(CXRB) to 1.5 x 10 -4 counts s _1 arcmin -2 found from 
the outer part of the image, and included 5% of the back- 
ground value as a systematic error, due to variation on the 
sky. We used XSPEC to convolve the surface brightness 
model through a spatial response matrix (constructed from 
the ROSAT PSF at 1 keV, for an azimuthally symmet- 
ric source centered on-axis) and to compare the convolved 
profile with the data. We find an acceptable fit in the ra- 
dial range 3'-43' (see Figure || and Table [l]) , with best fit 
parameters a x — 3.85 ±0.35 arcmin (= 310 ±30 h^ kpc), 
(3 = 0.68 ± 0.03 with X 2 = 74.1 for 87 degrees of freedom. 
Our values of a x and (3 are consistent with another study 
of the ROSAT PSPC data of A3571 (Vikhlinin et al. 1999) 
who also excluded the cooling flow area from the fit. In 
yet another study of ROSAT PSPC data of A3571 (Mohr 
et al. 1999) inconsistently smaller values were found for 
a x and (3. This is a consequence of including the data of 
the cooling flow into the profile fit in that work. 

If we assume that the intracluster gas is isothermal and 
spherically symmetric, the best-fit parameters a x and (3 
determine the shape of the gas density profile as: 



Pgas{r) = Pgas{0) 




(2) 



The observed temperature variation in A3571 from 7 keV 
to 4 keV with radius will introduce at most a 2% effect on 
the gas mass (e.g. Mohr et al. 1999), which is negligible 
compared to other components in our error budget. 

Wc obt ained the normalizatio n of the gas density pro- 
file (as in |Vikhlinin et al. 1999[) p gas (0) = 1.5 x 1O 14 M 



Mpc~ 3 , or 1.0 xl0~ ZD g cm - ' 3 , by equating the emission 
measure calculated from the above equation, with an ob- 
served value of 8.1 xlO 67 cm~ 3 inside a cylinder with 



r = 0.1 — 2 h§ Mpc radius, centered at the cluster bright- 
ness peak (r = 0.1 Mpc encompasses the cooling flow ex- 
cluded from all our analyses). 

3. TEMPERATURE DATA 

We used the temperature profile data presented in 
Markevitch et al. (1998), who combined three ASCA 
pointings to derive the emission weighted, cooling flow - 
corrected temperature kT = 6.9 ± 0.2 keV. Outside the 
cooling radius, the temperature values were measured in 
radial bins 2'-6'-13'-22'-35' (0.13-0.39-0.85-1.43-2.28 
Mpc). The central bin that is affected by a significant cool- 
ing flow component is not used in the analysis below. The 
temperature errors were determined by generating Monte 
- Carlo data sets which properly account for the statis- 
tical and systematic uncertainties (including those of the 
PSF, effective area and background). As in most nearby 
clusters, the ASCA data reveal a temperature decline with 
radius. The ROSAT PSPC data on A3571 in the 0.2-2 keV 
band were also analyzed by Irwin et al. (1999), who derive 
a temperature profile consistent with a constant up to 20'. 
However, those authors did not include the PSPC calibra- 
tion uncertainties that dominate the ROSAT temp erature 
errors for hot clusters such as A3571 (see e.g. Mark© 
vitch & Vikhlinin 1997a); inclusion of these uncertaintincs 



should make their results consistent with the ASCA pro- 
file. The ASCA temperature profile for A3571 is si milar to 
profiles of a large s ample of nearby ASCA clusters ( Marke 
vitch ct al. 1998), when scaled to physically meaningful 
units of the radii of fixed overdensity. Therefore it appears 
unlikely that the observed decline is due to an unknown 
instrumental effect. A more detailed discussion about the 
validity of the ASCA spatially resolved temperature data 
can be found in Nevalainen et al. (1999a). Hydrodynamic 
simulations predict a qualitatively sim ilar radial tcmpera- 
turc behavior in relaxed clust ers (e. g., lEvrard ct al 1996| 



Eke, Navarro, & Frenk 1997; pryan fc Norman 1997| ), al- 
though there are differences in detail between the simu- 
lations and observations as well as between the different 
simulation techniques (e.g., Frenk et al. 1999). 

4. VALIDITY OF THE HYDROSTATIC EQUILIBRIUM 

Quintana & de Souza (1993) report preliminary results 
of an optical study of the galaxies in A3571. They find 
a suggestion that the galaxy distribution in A3571 is ir- 
regular and forms several velocity subgroups, but they did 
not perform quantitative statistical analyses of the galaxy 
distribution, due to the small number of observed galaxies 
(see Figure 1 for the distribution of A3571 member galax- 
ies from the NASA Extragalactic Database). The central 
giant galaxy MCG05-33-002 has an extensive optical halo 
with dimensions of 0.2 x 0.6 h^Mpc, elongat ed alone 
the major axis of the core region of this galaxy ( Kemp & 
Meaburn 1991). The galaxy di stribution of A3571 is als o 
aligned in the same direction (Kemp & Meaburn 1991) 



As discussed by Quintana & de Souza (1993), the optical 
data suggests that the cD galaxy formed during the origi- 
nal collapse of the central part of the cluster and that the 
cluster may not yet be virialized. However, as Quintana 
& de Souza (1993) state, their galaxy distribution results 
are only tentative. Furthermore, galaxies are not the best 
measure of the relaxation since clusters form within in- 
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tersecting filaments in larger scale and superpositions can 
give the appearance of the asymmetries, substructure, and 
superposed groups. 

In X-rays, the ROSAT PSPC data of A3571 (Figure 1) 
show that the gas is azimuthally symmetric (except for a 
slight ellipticity, see Section ||) and that there is no sub- 
structure and no correlation between the galaxy and gas 
distributions at large radii. Furthermore, the ASC A gas 
temperature map of A3571 ( Markevitch et al. 1998 ) shows 
no asymmetric variation that would indicate dynamic ac- 
tivity. 

Neumann & Arnaud (1999) found evidence in their 
ROSAT cluster sample that cooling flows are a recurrent 
phenomena that may be turned off by mergers, in accor- 
dance with a hierarchical clustering scenario (Fabian et al. 
1994). Since A3571 has a considerable cooling flow (Peres 
et al. 1998), any merger must have been either not very 
strong or sufficiently in the past for the gas to reestablish 
equilibrium and a cooling flow. 

A3571 is a member of the Shapley supercluster (Ray- 
chaudhury et al. 1991) and therefore likely to have more 
frequent mergers and may not be typical of more isolated 
clusters. However, all the X-ray evidence consistently 
argues against any significant ongoing merger in A3571. 
Since the optical evidence does not contradict significantly 
the X-ray evidence for non-merger, we assume that the hy- 
drostatic equilibrium is valid in A3571. 

5. MASS FITTING 
5.1. Method 

For the details of the mass calculation, we refer to our 



similar analysis of cluster A401 ( Nevalainen et al. 1999a ). 
Briefly, we model the dark matter density with a constant 
core model 

r 2 \ -a/2 

Pdark OC \ 1 + — ) , (3) 



r 

l + — 

a d 



and with the central cusp profile: 



Pdark OC 



a,, 







)> 


--) 




ad J 



(4) 



We fix r\ = 1 in t he cusp models, as su ggested by numer- 
ical simulations (Navarro et al. 1997, hereafter NFW), 
but vary the other parameters. We solve the hydrostatic 
equilibrium equation 



M tot (<r) = 3.7OxlO 13 M - T(,) 



keV Mpc 



din Pgas dlnT 
d In r dhxr 
(5) 

(e.g. Sarazin 1988, using fx = 0.60), for temperature, in 
terms of dark matter and gas density profile parameters. 
We fix the gas density to that found from the ROSAT 
data above, calculate the 3 - dimensional temperature pro- 
file model corresponding to given dark matter parameters, 
project it on the ASCA annuli, compare these values to 
the observed temperatures and iteratively determine the 
dark matter distribution parameters. To propagate the 
errors of the temperature profile data to our mass values, 
we repeat the procedure for a large number of Monte - 
Carlo temperature profiles with added random errors. We 
reject unphysical models that give infinite temperatures 



at large radii, and those models that are convectively un- 
stable (that is, correspond to polytropic index > | in the 
radial range of the temperature data, outside the cool- 
ing flow region r = 3'-35'). From the distribution of the 
acceptable Monte - Carlo models, we determine the 1 a 
confidence intervals of the mass values as a function of ra- 
dius. We convert these values to 90% confidence values, 
assuming a Gaussian probability distribution. We cannot 
constrain all dark matter model parameters independently 
due to the limited accuracy of the temperature data. How- 
ever, the models with steeper dark matter density slopes 
(higher a) require larger dark matter core radii (higher a^) 
to produce similar shapes of temperature profile and due to 
this correlation the corresponding mass values vary within 
a relatively narrow range. We also propagate the estimate 
of the uncertainty of the local gas density gradient to the 
total mass values, as in Nevalainen et al. (1999a). 

In Figure |^ we show representative density profiles of 
forms (|j) and (^), and the corresponding model temper- 
ature profiles. Both functional forms give acceptable fits 
to the data and yield masses consistent within 90% con- 
fidence errors. Our final 90% confidence intervals of total 
mass, at each radius, include the 90% confidence intervals 
of both models, and the average of the two models is used 
as the best value. As can be seen in Figures p$|i and ||d, 
in the radial range 3' — 15', the Monte-Carlo densities are 
lower than the best fit values. This asymmetry is due to 
the fact that the polytropic index criterion effectively re- 
jects the most massive models, as was also found for A401 
(Nevalainen et al. 1999a). 

5.2. Results 

The final mass profile is shown in Fig. |]. The over- 
density, or the mean interior density in units of the criti- 
cal density, calculated from our best fit models, is 240 at 
r = 35', the lar gest radius covered by the ASCA data. 
Simulations (e.g. Evrard et al. 1996) suggest that within 
rsoo, where the overdensity is 500, hydrostatic equilibrium 
is valid. For A3571, 



r-500 - 25.9' = 1.7 h^ 1 



and the mass within this radius 



M tot (< r 5 oo) = x lO 14 ^ 1 M ( 



Mpc 



(6) 



(7) 



The mass values within several interesting radii are given 
in Table 1. 

At large radii our mass errors are quite large because 
they cover the values allowed by two different models, and 
include the uncertainty of 0. Therefore the isothermal 
mass is consistent with our results, but compared to our 
best values, the isothermal ones are greater by factors of 
1.1 and 1.3 at radii of r 50 o and 35' (see Figure Q). This dif- 
ference is a natural consequence of the real temperatures 
being lower than the average temperature at large radii, 
similarly with other clusters with measured temperature 
profiles (Markevitch fc Vikhlinin 1997b|, [Nevalainen et al 



1999a, Markevitch et al. 1999). However, at small radii 



[r = a x ), differently from the other above clusters, the 
isothermal mass in A3571 is about equal to the value ob- 
tained with the observed temperature profile, due to the 
nearly constant temperature up to 13' in A3571. 
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X-ray total mass estimate for A3571 



The deprojection method, with the isothermal assump- 
tion, gives a total mass value of 6.22 x lO 14 /^ 1 M in- 
side a radius of 0.91 Mpc ( |Ettori et al. 1997| ), 
whereas our value at that radius is significantly smaller, 
4.9±o g x lO 14 / ^ 1 M R . This behaviour is similar to that 



found for A401 (Nevalainen et al. 1999a). 

The frequently used mass - tempe rature scaling law ob- 
tained in cosmological simulations ( Evrard et al. 1996 ) 
predicts that A3571, which has kT = 6.9±0.2 kcV (Marke- 
vitch et al. 1998), will have r 500 = 2.1 ±0.14 M PC and 
M to t{< ?"50o) = 1.3±0.19 x lO 15 /^ 1 Mq, whereas our mea- 
sured values above are smaller by factors of 1.2 and 1.6. 
Similar behaviour has been found in several other hot clus- 
ters wi th temperature profiles measur ed with ASCA, i.e. 
A225 6 (iMarkevitch fc Vikhlinin 1997b|), A401 flNevalainciJ 



et al. |1999aD , A496 and A2199 flMarkevitch et al. 1999| ). 
Also temperature profiles of cooler groups NGC5044 and 
HCG62 and a galaxy NGC507, measured with ROSAT, 
give similar results (Nevalainen et al. 1999b). These com- 
parisons suggest that the above simulations produce too 
small temperature for a given mass. 

The virial theorem analysis of the galaxy velocity distri- 
bution in A3571 (Girardi et al. 1998) gives R vir = 4.18/i^ 1 
Mpc and M vir = 1.63t°;72 x 10 15 h^M Q , whereas our val- 
ues extrapolated to this radius are 1.2 ± 0.7 x W 15 h^M Q , 
consistent within the large errors. 



deed be universal. 

6.2. Gas mass fraction 

The dark matter density in the best fit model falls as r~ 4 
at large radii, whereas the gas density falls as r -2 . This 
causes the gas mass fraction to increase rapidly at large 
radii, to a value of 



/ 9 a S (<r 5 oo)=0.19t^^50 3/2 - (8) 
(see Fi gure 3). This value is consistent with those for 



A2256 (Markevitch & Vikhlinin 19971 


), A401 flNcvalainer 


|et al. 1999a|), A496 and A2199 (Mar 


kevitch et al. 1999) 



the gas and dark matter distributions are similar in dif- 
ferent clusters. Our value is also consistent with results 
for samples of clusters analyzed assuming isothermality: 
f gas (< r 500 ) = 0.168t°;^/i- 3/2 ( [Ettori fe Fabian 1999|) 



and f ga ,s(< f soo) — 0-212 ± 0.0 06 h R ^ 2 (for clusters with 
kT > 5 keV |Mohr et al. 19991). At larger radii the tem- 

Id\ 



perature profile analysis would probably yield still higher 
values compared to the isothermal ones. 

Following the method of White et al. (1993; see also 



Nevalainen et al. 1999a for its application to A401), us- 
ing the f gas value above, we compute an estimate for the 
cosmological density parameter Q m =< p > / p cr u as 



6. DISCUSSION 

6.1. NFW profile 

With i] = 1 and a = 3 the cusp model (Eq. 4) corre- 
sponds to the NFW "universal" mass profile. For A3571 
there is a significant detection of a cooling flow in the 
center and the polytropic 7 < | constraint is applica- 
ble only beyond the cooling flow region. At those radii 
the temperature gradient of the cusp model is not strong, 
and the model is convectively stable. This model also 
gives an acceptable fit to the A3571 data, and is consis- 
tent with our mass profile within the errors (see Figure 
|). Using the best fit NFW profile for A3571, we ob- 
tain the concentration parameter c = ?"2oo/ad = 5.3 and 
M200 — 1-3 x 10 15 h^Q M & . These values are consistent 
with the NFW simulations in SCDM and CDMA cosmo- 
logical models. 

In the hydrostatic equilibrium scheme, since the ob- 
served gas density and temperature profiles are similar in 
differe nt clusters, when sca led by their estimated virial 



s2 ), 



"f fga 



radii QVikhlinin et al. 1999| and [Markevitch et al. 199~8 



respectively), a similar total mass distribution is implied, 
and that is what we observe. The shapes of the total 
mass profiles at large radii in other hot cluster s with mea- 
sured ASCA tempe rature pr ofiles for A2256 (Markevitch! 
fc Vikhlinin 1997b| ), A2 029 (|Sarazin, Wise fc Markevitch! 
19981), A496 and A2199 (IMarkevitch ct al. 1999| ) and A401 
( Nevalainen ct al. 1999a] ) arc also consistent with the 



NFW model. NFW profile also describes well the mass 
profiles of cool groups NGC5044 and HCG62 and a galaxy 
NGC507 that are derived using ROSAT PSPC tempera- 
ture profiles (David et al. 1994, Ponman & Bertram 1993, 
Kim & Fabbiano 1995, respectively, see the discussion of 
NFW models for these objects in Nevalainen et al. 1999b). 
These consistencies suggest that the NFW profile may in- 



M, 



gal 



(9) 



at r5oo for A3571, where T is the local baryon diminution, 
for which we use a value 0.90 as suggested by simulations 
(Frenk et al. 1999). Due to the lack of a reliable esti- 
mate of the total galaxy mass in A3571 within rsoo we 
compute only the upper li mit for Q m using M„ n i = 0. We 
take n b = 0.076 ± 0.007 (|Burles ct al. 19981 ). Using a 



reasonable lower limit for the Hubble constant Hq > 60 
km s _1 Mpc -1 (e.g. Nevalainen & Roos 1998) we obtain 
fi m < 0.4, consistent with several independendent current 
fl m estimates (Freedman 1999, Roos & Harun-or-Rashid 
1999). 

7. CONCLUSIONS 

We have constrained the dark matter distribution in 
A3571, using accurate ASCA gas temperature data. The 
dark matter density in the best fit model scales as r -4 at 
large radii, and the NWF profile also provides a good de- 
scription of the dark matter density distribution. The total 
mass within r 500 (1.7 M Pc) is 7.8^2:2 x 1Ql4 h 50 M Q at 
90% confidence, or 1.1 times smaller than the isothermal 
value, or 1.6 times smaller than that predicted by the scal- 
ing law based on simulations (Evrard et al. 1996), which is 
qualitatively similar to the results for other clusters with 
accurate temperature profiles. The gas density profile in 
A3571, proportional to r~ 21 at large radii, is shallower 
than that of the dark matter. Hence the gas mass fraction 

increases with radius, with f gas {r5oo) = 0.19^0 03 h^^ 2 
(90 % errors) at r 50 o, consistent with results for A2256, 
A401, A496, and A2199. Assuming that this is a lower 
limit of the primordial baryonic fraction, we obtain Q m < 
0.4 at 90% confidence. However, f gas is still strongly in- 
creasing at rsoo , so that we obviously have not reached the 
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universal value of the baryon fraction, which would make 
f2 m even smaller. 
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6 



X-ray total mass estimate for A3571 



Table 1 
Mass of A3571 



Parameter Value 

a x a [arcmin] 3.85 ± 0.35 

a x a [h^ kpc] 310 ± 30 

13 a 0.68 ±0.03 

Pgas (0) [10- 26 h^ 2 g cm- 3 ] 1.0 

M(< 3.85' = a x ) [M Q ] 6.2j£J x 10 13 

M(< 15.2' = 1 Mpc) [M Q ] bAtil x 10 14 

M(< 25.9' = r 500 ) [M ] 7.8±H x 10 14 

M(< 35') [M ] 9.1±£! x 10 14 

f gas (< 25.9' = r 50 o) x h% 2 0.19±8;g| 



* using if = 50 /150 A:m s 1 Mpc 1 . All errors 
are 90% confidence errors. 

a From ROSAT PSPC data excluding the cen- 
tral r < 3' cooling flow region. 
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Fig. l.— ROSAT PSPC contour map of the surface brightness of A3571 in 0.73-2.04 keV energy range, smoothed by 
a Gaussian with a = V . The contour level values are 0.0005, 0.001, 0.002, 0.004, 0.008, 0.016 and 0.032 counts s _1 
arcmin" 2 . Galaxies from the NASA Extragalactic Database are plotted as filled circles. 
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X-ray total mass estimate for A3571 




Fig. 2. — ROSAT PSPC radial surface brightness profile in 0.73-2.04 keV range together with the PSF- convolved best-fit 
(3 model. The data with r < 3' are not included in the fit, due to central cooling flow. 



Nevalainen J., Markevitch M. and Forman W. 



9 



Core model Cusp model 
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Fig. 3. — In (a) and (b), crosses show projected ASCA temperatures with 1 a errors. Thin solid lines show a representative 
set of temperature models of the constant core form (a) and the cusp form (b), before projection, allowed by the convective 
stability constraint (see text). Thick solid line shows the best fit model (before projection). In the case of cusp model, the 
thick dotted line shows the values of this best model projected to the 2D ASCA bins, which are compared with the ASCA 
data. In (c) and (d), the corresponding density models are plotted, together with the gas density model. For comparison, 
values assuming isothermality (kT = 6.9 keV) are also shown. 
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X-ray total mass estimate for A3571 
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Fig. 4. — (a) shows the enclosed mass profile (solid line) with 90% confidence errors (shaded area), obtained combining the core and the cusp 
model fit results. Also the model assuming isothcrmality (dotted line), together with the gas mass profile (dash-dot) are shown. The best fit 
"universal" dark matter profile, suggested by NFW simulations, is plotted as dashed line. In (b), the gas mass fraction with errors is shown. 
Masses are evaluated using H = 50 km" 1 Mpc" 1 (total mass scales as H 1 and gas mass as 



